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We propose a QCD-like theory with the Zn^ symmetry. The flavor-dependent twisted boundary condition 
(TBC) is imposed on Nc degenerate flavor quarks in the SU(A^c) gauge theory. The QCD-like theory is useful to 
understand the mechanism of color confinement. Dynamics of the QCD-like theory is studied by imposing the 
TBC on the Polyakov-loop extended Nambu-Jona-Lasinio (PNJL) model. The TBC model is applied to two- 
and three-color cases. The Zn^ symmetry is preserved below some temperature Tc, but spontaneously broken 
above Tc. The color confinement below Tc preserves the flavor symmetry. Above Tc, the flavor symmetry is 
broken, but the breaking is suppressed by the entanglement between the Polyakov loop and the chiral condensate. 
Particularly at low temperature, dynamics of the TBC model is similar to that of the PNJL model with the 
standard fermion boundary condition, indicating that the Zjv^ symmetry is a good approximate concept in the 
latter model even if the current quark mass is small. The present prediction can be tested in future by lattice 
QCD, since the QCD-like theory has no sign problem. 

PACS numbers: 11.30.Rd, 12.40.-y 



I. INTRODUCTION. 

Understanding of the confinement mechanism is one of the 
most important subjects in hadron physics. According to Lat- 
tice QCD (LQCD), the system is in the confinement and chiral 
symmetry breaking phase at low temperature (T), but in the 
deconfinement and chiral symmetry restoration phase at high 
T. The confinement mechanism is, nevertheless, still unclear 
for several reasons. The main reason is that the exact sym- 
metry is not found for the deconfinement transition and hence 
the order parameter is unknown. In the limit of zero cuiTent 
quark mass, the chiral condensate is an exact order parame- 
ter for the chiral restoration. In the limit of infinite cuiTent 
quark mass, on the contrary, the Polyakov loop becomes an 
exact order parameter for the deconfinement transition, since 
the Zjv^ symmetry is exact there. For the real world in which 
u and d quarks have small current quark masses, the chiral 
condensate is considered to be a good order parameter, but it 
is not clear whether the Polykov loop is a good order param- 
eter. In this paper, we approach this problem by proposing a 
QCD-like theory with the Zat^ symmetry. 

We start with the SU(Afc) gauge theory with Nf degenerate 
flavor quarks. The partition function Z in Euclidean space- 
time is described by 

Z = J DqDqDAexpi-So] (1) 

with the action 

So^J d^xiY, qfil.D, + mf)qf + ^F^,\ (2) 
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where q/ is the quark field with flavor / and current quark 
mass nif, = + lA^ is the covariant derivative with the 
gauge field A^, g \s the gauge coupling and _F)j,y = — 
duA^, - i[A^,A^,] = F^^T" with the SU(7Vc) generator T". 
The temporal boundary condition for quark is 

qf{x,p = l/T) = -qf{x,0). (3) 

The Ztv^ transformation changes the ferimon boundary con- 
dition as HH 

/3) = - exp (-z2^/c/iV,)g/(x, 0) (4) 

for integer k, while the action 5*0 keeps the original form ^ 
since the 'Lm. symmetry is the center symmetry of the gauge 
symmetry (ill- The Zjv^ symmetry thus breaks down through 
the fermion boundary condition in QCD. 

Now we consider the SU(Af) gauge theory with N degen- 
erate flavor quarks, i.e. N = Nf = Nc, and assume the 
following twisted boundary conditions (TBC): 

qf{x,/3) = -exp{ief)qf{x,0) 

= -exp[t{0^+27T{f-l)/N)]qf{x,O) (5) 

for flavors / labeled by integers from 1 to A^; see FiglUfor 
the twisted angles 9f. Here 9i is an arbitrary real number in 
a range of < 6'i < 27r. The action 5*0 with the TBC is 
not QCD but a QCD-like theory. The QCD-like theory has 
the Zat^ symmetry, i.e. invariant under the transforma- 
tion. In fact, the Zat^ transformation changes / into f — k, 
but f — k can be relabeled by / since Sq is invariant under 
the relabeling. The QCD-like theory with the Zn^ symmetry 
is useful to understand the mechanism of color confinement. 
When the fermion field qf is transformed by 

qf ^ exp {—i9 fTT)qf (6) 



for Euclidean time r, the action Sq is changed into 




(7) 
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Fig. 1: Twisted factors e'^f (f = 1, 2, 
complex plane for the case of 9i = 0. 



, A'') on a unit circle in the 



with the imaginary quark number chemical potential /z/ = 
iTOf, while the TBC is transformed back to the standard one 
([3]l- The action Sq with the TBC is thus equivalent to the ac- 
tion S{9f) with the standard one ([3]). In the hmit of T = 0, 
the action S{6f) tends to 5*0 with OJ fixed. The QCD-like the- 
ory is thus identical with QCD at T = where the Polyakov 
loop # is zero. This indicates that in the QCD-like theory the 
Ziv^ symmetry is preserved up to some temperature and 
spontaneously broken above Tc- In the QCD-like theory, the 
flavor symmetry is explicitly broken by the flavor-dependent 
TBC. As shown later, the flavor-symmetry breaking is recov- 
ered below Tc. The breaking then becomes significant only 
above Tc- 

In general, the QCD partition function Z{T. 6) with finite 
imaginary chemical potential 6 has the Roberge-Weiss (RW) 
periodicity |1]: Z{T,e) = Z(T, 6* + 2'Kk/Nc) for any inte- 
ger k. The RW periodicity was confirmed by lattice QCD 
(LQCD) lHHH and the Holographic QCD lUl. The RW pe- 
riodicity means that Z{T, 9) is invariant under the extended 
Zjv^ transformation, i.e. the combination of the Zjv^ transfor- 
mation and the parameter transformation 6^9 + 2Trk/Nc. 
Actuafly, Z{T, 9) is transformed into Z{T, 9 - 2'Kk/Nc) by 
the Zjv^ transformation and Z{T, 9 — 2TTk/Nc) is transformed 
back to Z{T, 9) by the parameter transformation. The QCD 
partition function thus has the extended Zjv^ symmetry, and 
dynamics of QCD at imaginary chemical potential is governed 
by the symmetry. 

The extended Zjv^ symmetry is not an internal symmetry, 
since the transformation includes the shift of external param- 
eter 9. In the QCD-like theory, the shift of 9 is not nec- 
essary because of the TBC. Thus the QCD-like theory pos- 
sesses the Ztv^ symmetry as an internal symmetry, whereas 
QCD has the extended Zjv^ symmetry as an external sym- 
metry. The Polyakov-loop extended Nambu-Jona-Lasinio 
(PNJL) model |[2l|l4M5lll is a good model to understand QCD 
at finite imaginary chemical potential 9 and hence the QCD- 



like theory, since the PNJL model possesses the extended Z 
symmetry in the standard fermion boundary condition (O |2|]. 

In this paper, we propose a QCD-like theory with the 
Zjv^ symmetry. The theory is constructed by imposing the 
TBC on the SU(Afc) gauge theory with Nc degenerate fla- 
vor quarks. Dynamics of the QCD-like theory is studied con- 
cretely by imposing the TBC on the PNJL model. Two cases 
of Nc = Nf = 2 and 3 are mainly considered. In this pa- 
per, the PNJL model with the TBC is shortly called the TBC 
model, and the PNJL model with the standard boundary con- 
dition is named the standard-PNJL model. We first show that 
the symmetry is preserved below some temperature Tc, 
but spontaneously broken above Tc- The interplay between 
the Zat^ symmetry breaking and the flavor symmetry break- 
ing is investigated. Comparing the deconfinement transition 
in the TBC model with that in the standard-PNJL model, we 
show that the Z^r^ symmetry is a good approximate concept in 
the latter model, even if the cuiTent quark mass is small. The 
present prediction can be checked by LQCD in future, since 
LQCD with the TBC is free from the sign problem. 



This paper is organized as follows. The case of Nc 



N 



f 



2 is investigated in Sec. HJand that of Nc = Nf = 3 is in Sec. 
Uni Two interesting extensions of the TBC model are shown 
in Sec. IIIICI Section UVj is devoted to summary. 



II. CASE OF iVc = 2 

A. Formalism 

The two-color and two-flavor PNJL Lagrangian ll34ll in Eu- 
clidean spacetime is 



(1 - a)Gs + {qfil5Taqff] 

f 1=0 



det (g,(l+75)gj) + det(gi(l 



+'iaG: 
+UmAmA]* ,T) 



15) qj) 



(8) 



with — d„ + iA„ — dy + iSuAA^^a-Y for the gauge 
field AJ^, where the {fa) for a = 1,2,3 are the Pauli 
matrices in flavor (color) space and tq is the unit matrix in 
flavor space. In the NJL sector, (1 — a)Gs denotes cou- 
pling constants of scalar- and pseudoscalar-type four-quark 
interactions, whereas aGs is that of the Kobayashi-Maskawa- 
't Hooft determinant interaction [52, 531 . Here a can vary 
from to 1/2 for positive Gg- The J7a(1) anomaly vanishes 
when a = 0. The Polyakov potential U, defined in ( fTTT l. is 
a function of the Polyakov loop (p and its Hermitian conju- 
gate <P*. The parameter TO/ (^/) stands for the cuiTent quark 
mass (the chemical potential) for each flavor. Here we set 
mo = TUu ^ md- 

In the PNJL model, the gauge field A^ is treated as a ho- 
mogeneous and static background field [l^S]- In the case 
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of Nc = 2, the Polyakov-loop <P and its conjugate ^* are de- 
termined in Euclidean spacetime by 



^=^tr,(L), <?* = itr,(Z), 



(9) 



where L = cxp{iA4/T) with A4 = {Aq. In the 
Polyakov-gauge, A4 is diagonal in color space, i.e., A4/T = 
diag(0i, (P2) for the (pi satisfying (j>i + (j>2 — 0. This leads to 



<p = 










e-*^i) = cos((^i), 


= 










+ e*^i) = cos(0i) 



(10) 



indicating that ^ is real. For the Polyakov-loop potential U, 

we use 



U = -6T[24e-''/'^^>2 + log (l - $^)] 



(11) 



proposed in Ref. fH, where a = 858.1 MeV and b^/^ = 
210.5 MeV. The Polyakov potential yields the second-order 
deconfinement phase transition at Tc — 270 MeV in the pure 
gauge theory. 

Now we consider the imaginary chemical potential /i/ = 
i9fT, where the twisted angles 6'/ are real. Making the 
mean-field approximation (MFA) and the path integral over 
the quark fields in the PNJL partition function .^pnjl. one 
can obtain the thermodynamic potential (per unit volume) as 



-Tln(ZpNjL)/V^ 

f=u,dc=l,2 



(27r)3 



Ef 



^In [1 
iln [1 



+ U{a,ao)+ll{^,T), 
where Ef{p) = Ef{p) ± fif foiEf{p) = 



(12) 



Mf 



Mu = 7no - 2Gs(cr + Cflo), 
Md = Too - 2Gs(cr - Cflo), 



(13) 
(14) 
(15) 



C = 1 — 2a, a = {uu + dd) and ao = {uu — dd). Here 
only the flavor-diagonal scalar condensates are taken. On the 
right-hand side of ( fT2] i only the first term is regularized by 
the three-dimensional momentum cutoff A lfl6l [itIi . since it 
diverges. 

The variables, X = ct, ao), are determined by the 

stationary conditions 



dSl/dX = 0. 



Solutions X{T, Of) of the conditions do not give a global min- 
imum of [2 necessarily, when the solutions are inserted back 
to (fT2l i. There is a possibility that they yield a local minimum 
or even a maximum. We have then checked that the solutions 
yield a global minimum 

Following Ref. ll34ll . we take the parameter set of toq — 
5.4 MeV, A = 657 MeV and Gs = 7.23 GeV^ that yield 
— (Sw)^/3 = 218 MeV, the pion decay constant /^r — 
75.4 MeV and the pion mass to^t — 140 MeV at vacuum. 

Taking the summation over color indices in fT% leads to 



[2 



d^ 



_2 y / ^ 2Ef 

+ ^ln [1 + G2,i(p)e^'^ + C2,2(p)e'''^ 



In [1 + G2,i(p)e-^'^ +G2,2(p)e- 



where 



G2,i(p) = 2<^'e-''^^ 
G2.2(p) = 



(17) 



(18) 



It is found from (O that G2.1 = and G2,2 7^ when ^ = 0. 
The configuration means that two colored quarks are statisti- 
cally in the same state. The configuration is thus realized as a 
result of the color confinement (^ = 0). In other words, the 
color confinement can be defined by the configuration. 
Making the Z2 transformation 



in ST7\ . one can find that i7 has the RW periodicity. 



(19) 



(20) 



(16) 



Namely, f2 has the extended Z2 symmetry. The TBC corre- 
sponds to setting 9^ = Oi and 9d ~ Oi+tt in Q. The Z2 sym- 
metry with odd k changes 6d) to {Od, Ou), but {9d, 9^) re- 
turns to (^u, 6d) by the relabeling of flavors. The TBC model 
thus has the Z2 symmetry as an internal symmetry in addition 
to the extended Z2 symmetry as an external symmetry. 

In the color-confinement phase defined by = 0, the ther- 
modynamic potential Q has only the configuration of G2,i = 
and G2,2 7^ 0, as mentioned above. Components includ- 
ing G2.2 in ([TtI i have flavor dependence only through factors 
g±2ie/^ but the factors has no flavor dependence because of 
{Ou, Od) = {01, Oi + tt). Noting that the M / are determined 
by the stationary condition (fTSI l from the flavor-independent 
J7, one can see that the flavor symmetry is recovered by the 
color confinement. 



B. Numerical results 

Let us start with the standard-PNJL model, i.e., the PNJL 
model with no chemical potential. In this case, the Polykov 
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loop ^ is an approximate order parameter of the color confine- 
ment, while the chiral condensate a is an approximate order 
parameter of the chiral transition. The flavor symmetry break- 
ing is described by the isovector condensate oq. We mainly 
consider the ?7a(1) symmetric case by taking a = 0. 

Figure |2] shows T dependence of and cr calculated with 
the standard-PNJL model, where a is normalized by cto = 
a{T — 0, fif =0). Both cr and are finite for any T, since 
there is no exact chiral and Z2 symmetry. As T increases, a 
decreases gradually, while <1> increases smoothly. The chiral 
and deconfinement transitions are thus crossover. Here uq is 
zero at any T, since the flavor symmetry is not broken. 




Fig. 3: Twisted factors e f on a unit circle in the complex plane for 
the case of Si = (left) and 6*1 = -f (right). 
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Fig. 2: 

(0,0), 



0.1 0.15 0.2 0.25 
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T dependence of a and <l> in the case of a = and {8u,0d) 
Here a is normalized by ao = a{T = 0, /i/ = 0). 



Now we consider the TBC model with the Z2 symmetry. 
The Polyakov loop <P is an exact order parameter of the color 
confinement. When <P ^ 0, there are two Z2 vacua. The 
vacuum with positive ^ is taken in this paper 

First we analyze the case 



9rf) = (-V2,7r/2) 



(21) 



corresponding to 9i — — 7r/2 in the TBC of Q; see the right 
panel of Fig. |3]for the twisted angles. In this case, the flavor 
symmetry is not broken by the TBC, because 



(22) 



where the first and second equalities are obtained by the 
charge-conjugation and (l2Tl i. respectively. 

Figure |4] shows a and as a function of T; note that oq is 
zero for any T because of the flavor symmetry. The Polyakov 
loop ^ is zero up to T = w 260 MeV, but finite above 
Tc- The Z2 symmetry is thus preserved exactly below Tc, but 
spontaneously broken above Tc- The deconfinement phase 
transition is second-order, since <P has no jump at T = Tc- 
Meanwhile, the chiral transition is crossover There is no qual- 
itative difference between the standard-PNJL model and the 
TBC model with 9i = —tt/2 for the deconfinement and chiral 
transitions, although the order of the deconfinement transition 
becomes stronger by the exact Z2 symmetry. 

In Fig. |5] the color state factors 6*2,1 (p = 0) and C2.2(p — 
0) are drawn as a function of T. The one-quark state C2.1 (p = 
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T(GeV) 

Fig. 4: T dependence of a and <l> in the case of a = and {6u,0d) 
(— 7r/2, 7r/2). a is normalized by cro. 



0) vanishes below Tc because of ^ = 0. Above Tc, on the con- 
trary, the system is dominated by the one-color state, although 
the two-quark state C2,2 remains there. 



U 




0.2 0.3 

T(GeV) 



0.4 



Fig. 5: T dependence of the color state factors C2,i(p = 0) (solid 
line) and C2,2(p = 0) (dashed line) in the case of a = and 
{eu,ei) = (-7r/2,7r/2). 



The delay of the chiral restoration at higher T can be under- 
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stood as follows. Taking the flavor summation in ( fT2] l leads to 



n 



(27r)3 



In [l + F2,i(p)e-^=+F2,2(p)e-2'*=] 



+ Uia, a„)+U{$,T), 



(23) 



where 



F2AP) = e'''--^^' 
-F2,2(p) = e"^"" ■ ""'e 



i(eu+9d)„-;3(£;„+_Bd) 



(24) 



Since = 6*1 and = ^'i + tt, Eq. (l24b is reduced to 



^2,1 (p) 

^^2,2(P) 



(2:2,16 



22,26 



_g2i9ig-/3(£;„+Bd) 



(25) 



where Z2,i = 1 and Z2.2 = — 1 are elements of the Z2 group. 
In the case of {Ou,Od) = (— 7r/2, tt/2), the flavor symmetry is 
not broken, so that E = Eu = Ed- In this situation, i^2,i and 
F2 2 are further reduced to 



^^2,1(P) 
^^2,2(P) 



(22,1 + 22,2) e^''^ 

-2I3E 







(26) 



The thermodynamic system thus has no ^2,1 but finite ^2.2. 
This means that u- and d-quarks are statistically in the same 
state. The chiral condensate a has weak T dependence, since 
the two-quark state factor F2.2 is strongly suppressed by the 
factor cxp( — 2f5E). Eventually, the chiral restoration becomes 
much slower in the TBC model. This slow restoration is true 
also for the case of iVc = 3 and A^/ = 2 |'39'], although the Z3 
symmetry is not exact in the case. 
Next we analyze the case 



{Ou^ed) = (0,^) 



(27) 



corresponding to 0i = in the TBC of (|5]i; see the left panel 
of Fig. |3] for the twisted angles. Figure |6] presents T depen- 
dence of flo and (p. In this case, the flavor symmetry is ex- 
plicitly broken by the TBC. The second-order deconfinement 
phase transition occurs at T = Tc ~ 235MeV. Below Tc, ao 
and (p are zero, indicating that the flavor symmetry is restored 
by the color confinement. Above Tc, both oq and become 
finite, indicating that the flavor and Z2 symmetries break si- 
multaneously. At high T where the flavor symmetry breaking 
is strong, a is getting large with respect to increasing T. This 
behavior is quite different from the corresponding behavior of 
a in the standard-PNJL model. 

Figure |7] shows T dependence of the constituent quark 
masses AI^ and M„. The quark masses are degenerate be- 
low Tc, but above Tc d-quark becomes heavier while w-quark 
does lighter The mass splitting is a consequence of the flavor 
symmetry breaking. 
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T(GeV) 

Fig. 6: T dependence of cr (solid line), ao (dashed line) and ^ 
(dotted line) in the case of a = and (Ou, 0d) ~ (0, tt). a and ao 
are normalized by ao- Note that a < and ao > 0. 




0.15 0.2 0.25 0.3 0.35 

T(GeV) 

Fig. 7: T dependence of the constituent quark masses in the case of 
a = and 6d) = (0, tt). The solid (dashed) line represents u 
(d) quark mass. 



In Fig. [8] the color state factors C2,i(p — 0) and C2.2(p = 
0) are plotted for w-quark as a function of T. Below Tc, only 
the two-quark state C2.2 remains. Above Tc, the system is 
dominated by the one-quark state C2,i. 

Figure |9] shows T dependence of ag and <P in the case of 
a = 0.2. The T dependence is similar to that in the case of 
a = 0, although Tc w 265 MeV in the former and 235 MeV in 
the latter. Comparing Fig. |9]with Fig. |6l one can see explicitly 
that Tc becomes larger as a increases. The C/a(1) anomaly 
thus delays the spontaneous breaking of the Z2 symmetry and 
hence that of the flavor-symmetry breaking. 



III. CASE 0¥ Nc = 3 

A. Formalism 

The present formulation for Nc = Nf = 3 is pararell to that 
for Nc = Nf = 2 shown in SecUlA] The PNJL Lagrangian 
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Fig. 8: T dependence of the color state factors C2,i(p = 0) (solid 
line) and 6*2,2 (p = 0) (dashed line) for u quark in the case of a = 
and (6»„,6ld) = (0,7r). 
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T(GeV) 

Fig. 9: T dependence of o (solid line), ao (dashed line) and <!> 
(dotted line) in the case of a = 0.2 and {Ou,0d) = (0, tt). a and ao 
are normalized by ao. Note that a < and ao > 0. 



with Nc = Nf = Sis 

£ = ^ qfiluD^ - Hf-fi + mf)qf 



f 1=0 

Gr,\detqi{l +75)9^ + detqi(l - 75)9^- 

L ij ij 



-U{<l>[Al<P[A],T), 



(28) 



where — du + iA^ ^ dy + i^y^^A^^a^aj"^ with the 
Gell-Mann matrices \a in color space. In the interaction part, 
Aa (a ^ 0) and Aq are the Gell-Mann matrices and the unit 
matrix in flavor space, respectively, and Gs and Gd are cou- 
pling constants of the scalar-type four-quark and the KMT de- 
terminant interaction ll52il53ll . respectively, in which the deter- 
minant runs in flavor space. The KMT determinant interaction 
breaks the U pSX) symmetry explicitly. 

The Polyakov-loop <P and its conjugate ^* are determined 



by 



7trc(i), 



(29) 



where L = exp(ij44/T) with A^jT = diag(0r, ^h)- 
Noting that (j)^ + 4>g + (j)b = 0, one can obtain 



J<l>g 



') 
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We take the Polyakov potential of Ref. 112011 : 



(30) 



h{T) ln(l - 6<P^'* + 4(<?>^ + <Z>*^) - 3(<P^>*)2) 



a(r) = ao + ai (^) + 02 (^) , biT) = 63 



(31) 



(32) 



Parameters of are determined to reproduce LQCD data at 
finite T in the pure gauge limit. The parameters except Tq 
are summarized in Table U The Polyakov potential yields the 
first-order deconfinement phase transition at T = To in the 
pure gauge theory Il54i |55[1 . The original value of Tq is 270 
MeV determined from the pure gauge LQCD data, but the 
PNJL model with this value of To yields a larger value of the 
pseudocritical temperature Tc at zero chemical potential than 
Tc « 160 MeV predicted by full LQCD |56-58J. We then 
rescale To to 195 MeV to reproduce Tc = 160 MeV Q. 



ao 


ai 


a2 


b3 


3.51 


-2.47 


15.2 


-1.75 



TABLE I: Summary of the parameter set in the Polyakov-potential 
sector determined in Ref. 12(11 . All parameters are dimensionless. 



Now we consider the flavor-dependent imaginary chemical 
potential i^f — iSfT. The thermodynamic potential (per vol- 
ume) based on the mean-field approximation is 1,421 



(2^ 



Ef 



-2 E E 

f—u^d.s c—r,g,b 

P 

U{au,ad,as)+U{<P,<P*,T) 



(33) 
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with af = (qfqf) and Ej = y + A//^ for / — u,d,s, 
where the three-dimensional cutoff is taken for the momentum 
integration in the vacuum term I.4Z1 . The dynamical quark 
masses M f are defined by 

M/ = m/ - 4G'sCT/ + 2GDcr/'cr/" (34) 

for f f and / ^ /" and /" ^ /"'. The mesonic potential 
U{cru,crd,crs) are obtained by 

U{au,ad,<Js) ^ ^ 2Gsa) - AG-DUufTdCTs- (35) 

For the 2+1 flavor system with = = to/, the PNJL 
model has five parameters (Gs, Gd, m-z, "^s, A typical set 
is obtained in Ref. |59]. The parameter set is fitted to empiri- 
cal values of 77'-meson mass and vr-meson mass and 7r-meson 
decay constant at vacuum. In the present paper, we set to^ 
to mi in the parameter set of Ref. |59], since we consider the 
three degenerate flavor system with toq = to; = m^. The 
parameter set is shown in Table HIl 



mo(MeV) 


A{MeV) 


GsA^ 


Gd/1' 


5.5 


602.3 


1.835 


12.36 



TABLE II: Summary of the parameter set in the NIL sector. All the 
parameters except mo are the same as in Ref. jsgll . 

Taking the color summation in ( |33] ) leads to 

j—u,d,s 

+ iln[l + G3,i(p)e''^' 

+ iln[l + G*i(p)e-''/e-^^^ 

+ U{au,ad,(T,)+U{$,$*,T), (36) 

where 

G3.i(p) = 3<Z'e-^^^, 
C3.2(p) - 3<Z'*e-2^^/, 
CzAp) = e^'^'"'- (37) 
One can find that has the RW periodicity, 

Q{ef) = n{df + 2fc7r/3), (38) 
making the Z3 transformation, 

in n. In the case of (6l„,6ld, 61,) = (6li, 6*1 +27r/3, 6*1 +47r/3), 
f2 is invariant under the Z3 transformation, indicating that f2 



possesses the Z3 symmety. When <P vanishes, the color con- 
finement (G3.1 — G3^2 — 0) occurs and thereby i7 has the 
flavor symmetry (Eu = Ed = Eg) since the factors e^^^^f 
have no flavor dependence in ( |36] |. The flavor symmetry is 
thus preserved by the color confinement (^ — 0) also for the 
case of Nc = 3. 



B. Numerical results 

First we consider the standard fermion boundary condition 
by setting 0^ ^ Od = Og ~ 0. In this case, the cr/ are de- 
generate and hence a = (a„ + ct^ + o's)/3 = cr/. Figure [TOl 
shows T dependence of a and (p. Both the chiral restoration 
and the deconfinement transition are crossover, although the 
former transition is a bit slower than the latter [Isot] . 




0.1 0.15 0.2 0.25 0.3 



T(GeV) 

Fig. 10: T dependence of a and $ at 9u = 9d = 9s = 0. a is 
normalized by ao- 

Next we consider the TBC model by taking two cases of 
(On, Od, 9s) = (0, 27r/3, 47r/3) and (-tt, ~tt/3, 7r/3) that cor- 
respond to the left and right panels in Fig. (TT] respectively. 
The charge conjugation yields the relation 

n{9^,dd,9s) = f2{-9^,-9d,-9,) = f2i9u,9,,0d) (40) 

for the two cases. Thus s-quark is symmetric with d-quark in 
these cases. Because of the Z3 symmetry, there are three Z3 
vaccua when <P Q. We then take the solution in which a 
phase (p of (p lies in a range of — 7r/3 < (j> < Tr/3. In the 
solution, is found to be real. 

Figure [T2I shows T dependence of several physical quanti- 
ties in the case of {9^, 9d, 9s) — (— tt, — 7r/3, 7r/3). The order 
parameters <P, a and aq = (t„ — ad — (Ju — ^s are plotted in 
panel (a). The first-order deconfinement transition takes place 
at T = Tc « 195 MeV. Below T^, ao and ^> are zero. The 
flavor symmetry is thus preserved by the color confinement. 
Above Tc, ao and ^ become finite, indicating that the flavor 
and Z3 symmetries break simultaneously. 

For the constituent quark masses Mf shown in panel (b), 
all the Mf are degenerate below Tc- Above Tc, Mu becomes 
heavier whereas two of the three, Md and Ms, are degenerate 
and becomes lighter The increase of Af„ makes the chiral 
restoration slower In panel (c), the absolute values of the 
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U 




Fig. 11: Twisted factors f on a unit circle in the complex plane 
for the case of = (left) and 6i = — tt (right). 



color-state factors C3.1, C3.2 and C3.3 are plotted at p = 0. 
Below Tc, C3 3 is small but finite, whereas C3 1 — C3 2 = 
0. Above Tc, the system is dominated by the one-quark state 

Here the case of {du,0d,0s) — (0, 27r/3, 47r/3) is consid- 
ered briefly. As shown in Fig. [13] below Tc « 195MeV phys- 
ical quantities have the same properties as those in the previ- 
ous case. The difference between the two cases appears above 
Tc. Particularly for Mf, it is found that = Ms > Mu 
in the present case, while Md — Ms < in the previous 
case. Thus both d- and s-quarks becomes heavier as T in- 
creases from Tc- This property makes the chiral restoration 
even slower in the present case. 



C. Two extensions of the TBC model 

In this subsection, we extend the TBC model in two direc- 
tions. 

As the first extension, we use the entanglement PNJL (EP- 
NJL) model jH instead of the PNJL model. A possi- 
ble origin of the four-quark vertex Gs is a gluon exchange 
between quarks and its higher-order diagrams. If the gluon 
field Ai, has a vacuum expectation value (Ao), is cou- 
pled to (Aq) and hence to (p through L li6Qrl . This effect al- 
lows Gs to depend on ^, namely Gs = Gs(<?) IS- It is 
expected that <!> dependence of Gs (^) will be determined in 
future by accurate methods such as the exact renormalization 
group method ll60l - l62ll . In this paper, however, we simply as- 
sume the following Gs (^) by respecting the chiral symmetry, 
the charge-conjugation symmetry li37J and the extended Z3 
symmetry |2|]: 



Gsi-P) Gs[l - ai<P'P* - a2(<? + )]■ 



(41) 



The PNJL model with the entanglement vertex W\\ is called 
the EPNJL model lHHIl^. In principle, Gd can depend on ^, 
too. However, (^-dependence of Gd yields qualitatively the 
same effect on the phase diagram as that of Ggjp. We then 
neglect (^-dependence of Gd, following Ref. |50]. 

The parameters ai and 02 in (1411 1 are so determined as to 
reproduce two results of LQCD at finite T. The first is a result 
of 2+1 flavor LQCD at ^ = |63] that the chiral transition 
is crossover at the physical point. The second is a result of 
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Fig. 12: T dependence of (a) order parameters ct, ao, (b) con- 
stituent quark masses Af/, (c) color-state factors Ca.i, Cz,2, 6*3,3 at 
p = in the case of Qd, 9a) ~ (— tt, — tt/S, vr/3). Here a and 
ao are normalized by ao- Note that Aid ~ Ms, a < and ao < 0. 



degenerate three-flavor LQCD at 6* = tt llOtl that the order 
of the RW endpoint is first-order for small and large quark 
masses but second-order for intermediate quark masses. The 
parameter set (cti, 02) satisfying these conditions is located 
in the triangle region 115011 

{-1.5ai + 0.3 <a2< -0.86ai + 0.32, as > 0}. (42) 

As a typical example, we take ai — 0.25 and a2 = 0.1, 
foflowing Ref. iH and rescale To to 150MeV |50]. 

Figure [14] shows T dependence of a, gq and <P calculated 
with the EPNJL model for (a) {Ou, 9d, ds) = (0, 0, 0) and (b) 
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T(GeV) 




0.1 0.15 0.2 0.25 0.3 

T(GeV) 

Fig. 13: T dependence of (a) order parameters a, ao, (b) con- 
stituent quark masses A//, (c) color-state factors Ca,!, C3,2, Ca.s at 
p = in the case of (6u, 9d, 9s) ~ (0, 2n/3, 47r/3). Here a and ao 
are normalized by ao- Note that Aid ~ Ms, a < and ao > 0. 



{0u,0d,ds) = (0, 27r/3, 47r/3). In panel (a), the chiral restora- 
tion and the deconfinement transition are first-order, because 
of the small current quark mass (5.5MeV) and the strong cor- 
relation between a-f and ^ |50]. In panel (b), the TBC model 
with the entanglement vertex yields similar T dependence to 
the EPNJL model with the standard quark boundary condition 
for the chiral restoration and the deconfinement transition, al- 
though the flavor symmetry is broken above Tc- 

As the second extension of the TBC model with Nf = Nc, 
one can consider the TBC model with Nf = INc for any posi- 
tive integer /. It is obvious that the TBC model with Nf = INc 
has the Zjv^ symmetry, if the twisted angles Of sae properly 
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(b) ao 
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Fig. 14: T dependence of order parameters a, ao and ^ calcu- 
lated with the EPNJL model for (a) (S^, ^d, Qs) = (0, 0, 0) and (b) 
(Su, Qd, ds) — (0, 2tt/3, 47r/3). Here a and ao are normalized by 
(To- Note that ao = in panel (a) and ao > in panel (b), while 
(T < in both panels. 



ordered; for example, 

9f=9, + 27rif-l)/Nf, 

or 

9f = ei+2Tr{f-l)/N,, 
for/ = l,2,.- - ,7V/. 



(43) 



(44) 




Fig. 15: Twisted factors f on a unit circle in the complex plane in 
the case of A'^c ~ 3, Nf — 6 and 9i = n/6. In the left and right 
panels, the e'^f are obtained by ( 143 1 and ( 144b . respectively. 
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Let us consider the case of Nc = 3, Nf = 6 and 9i = tt/6. 
In Fig. [15] the left and right panels show the twisted angles 
defined by (03]) and (l44l i. respectively. Here we take the right- 
panel case as an example. The thermodynamic potential i? has 
the same form as ( |36] |. except the flavor summation is taken 
from / = 1 to 6. It is straightforward to show that J? has the 
RW periodicity and the Z3 symmetry. 

We take the same parameter set as in the case of Nf — 
Nc — 3, except Gg is taken as Gs — Gs,3 — Gd,3O'/(0)/2 = 
2.226 GeV^, where Gs.3 and Gd,3 mean Gg and Gd in the 
case of Nc = Nf = 3, respectively, and cr/(0) stands for af 
at T and 6*/ = in the case of Nc = Nf = 3. We 
keep the Polyakov potential U of ( |32] |. but neglect the KMT 
determinant interaction just for simplicity. 

Figure |6] presents T dependence of ct/ and <P for the right- 
panel case of Fig. [15] Below Tc « 190 MeV, the flavor sym- 
metry ((Ti = CT2 = • ■ • = cTe) is preserved by the color con- 
finement (<1> = 0). Above Tc, the flavor and Z3 symmetries 
break simultaneously. The flavor symmetry breaking is par- 
tial because / = 1 is symmetric with f = A, f — 3 with 
/ = 6, and f — 2 with / = 5. As a consequence of this 
property, the ct/ spht into three doublets. 
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1 ^.^ 


- 1,4 
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<I> 
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0.1 0.15 0.2 0.25 0.3 

T(GeV) 

Fig. 16: T dependence of a/ and <P in the right-panel case of Fig.[T5] 
The solid, dashed and dot-dashed lines represent tJi = 0-4, 0-3 = as 
and a2 ~ os, respectively, whereas the dotted line corresponds to 



IV. SUMMARY 

We have proposed a QCD-like theory with the sym- 
metry. The QCD-hke theory is constructed by imposing 



the flavor-dependent twisted boundary condition Q on the 
SU(Afc) gauge theory with Nc degenerate flavor quarks. Dy- 
namics of the QCD-like theory has been studied by impos- 
ing the TBC on the PNJL model. The TEC model has the 
Ztv^ symmetry and hence the Polyakov loop becomes an ex- 
act order parameter of the deconfinement transition. The TBC 
model is a good model to investigate the mechanism of color 
confinement. 

For both cases of Nf = Nc — 2 and 3, the Polyakov 
loop is zero up to some temperature Tc, but becomes finite 
above Tc- The symmetry is thus preserved below Tc, but 
spontaneously broken above Tc- Below Tc, the color confine- 
ment preserves the flavor symmetry. Above Tc, meanwhile, 
the flavor symmetry is broken explicitly by the TBC. The 
flavor-symmetry breaking makes the chiral restoration slower, 
but the entanglement interaction between a and <P makes the 
restoration faster The entanglement interaction thus sup- 
presses the flavor symmetry breaking. In the standard-PNJL 
model with degenerate flavor quarks, becomes finite but 
small at T lower than the pseudo-critical temperature, while 
the flavor symmetry is preserved. Dynamics of the TBC 
model is thus similar to that of the standard-PNJL model be- 
low Tc- The similarity is relatively worse above Tc, but it is 
improved by the entanglement interaction. One can then ex- 
pect that QCD with the approximate Z3 symmetry is similar 
to the QCD-like theory with the Z3 symmetry and hence that 
the Z3 symmetry is a good approximate concept in QCD, even 
if the current quark mass is small. 

The model prediction mentioned above can be tested with 
LQCD, since LQCD with the TBC has no sign problem. The 
QCD-like theory is useful to understand the mechanism of 
color confinement, since the Zat^ symmetry is exact. For ex- 
ample, it is quite interesting to see T dependence of the po- 
tential between q and q in LQCD with the TBC. 
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